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Abstract 

We consider a quartic interaction of the Higgs and Goldberger-Wise fields, which con¬ 
nects the mechanism of the extra dimension size stabilization with spontaneous symmetry 
breaking on our brane and gives rise to a coupling of the Higgs field to the radion and its 
KK tower. We estimate a possible influence of this coupling on the Higgs-radion mixing 
and study restrictions on model parameters from the LHC data. 


1 Introduction 

In the present paper we consider an extension of the Standard Model (SM) based on the 
Randall-Sundrum model with two branes stabilized by a bulk scalar held [D [21 [31 H, which is 
necessary for the model to be phenomenologically acceptable. A characteristic feature of this 
extension is the presence of a massive scalar radion held together with its Kaluza-Klein (KK) 
tower. These helds have the same quantum numbers as the neutral Higgs held. Thus, the 
radion held and its excitations can mix with the Higgs held, if they are coupled. 

Originally, a Higgs-radion coupling in the unstabilized Randall-Sundrum model arising due 
to a Higgs-curvature term on the brane was put forward in [5]. Then, such a coupling and the 
resulting Higgs-radion mixing in the case of the stabilized model were discussed in paper [6] 
without taking into account the KK tower of higher scalar excitations. The phenomenology of 
the Higgs-radion mixing originating from the Higgs-curvature term was also considered in view 
of the discovery of the Higgs-like boson at the LHC [7]; various assumptions about the masses 
and the mixings of the scalar states have been analyzed in papers [8]-[T6]. In particular, it was 
shown that the light radion-dominated state with mass below or above the observed 125-GeV 
boson is still not completely excluded by all the electroweak precision constraints and the LHC 
data. 

Here, we discuss a different mechanism of Higgs-radion mixing immanent in stabilized brane- 
world models, where a Higgs-radion coupling naturally arises due to spontaneous symmetry 
breaking on the brane involving the stabilizing scalar held. This approach takes into account 
the inhuence of the KK tower of higher scalar excitations on the parameters of the Higgs-radion 
mixing, which turns out to be of importance. 

In principle, the most general mechanism of Higgs-radion mixing can include both possible 
Higgs-radion couplings. However, after the spontaneous symmetry breaking on the brane the 
Higgs-curvature term also gives rise to a brane-localized curvature term that affects the mass 
spectrum and the couplings to matter of the graviton KK modes mi; this effect should be taken 
into account when examining the four-dimensional effective theory. For this reason, here we 
restrict ourselves to the new mechanism of Higgs-radion mixing, which has the advantage that 
it modihes only the scalar sector of the model and leaves intact the masses and the coupling 
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constants of the graviton KK excitations. This also allows one to isolate the effects dne only 
to the new mechanism of the mixing. 


2 Higgs-radion interaction 

A stabilized brane-world model in hve-dimensional space-time E = M 4 IZ 2 with coordinates 
{x^} = M = 0,1, 2, 3,4, = 0,1, 2, 3, with the coordinate = y, —L < y < L 

parameterizing the hfth dimension, is dehned by the action 

S = Sg + S^+SM, (1) 


where Sg and are given by 


= 


-2M^ j d^x j dyRy/^, 


S. 


<I)+SM 


'-L 


d X I dy{-g dM(j)dN(j) -V{(j))] ^/g 


^/^Vl{(j))d'^X+ y/^{-V2{(j)) + LsM)d‘^x. 


ly=0 


ly=L 


( 2 ) 

(3) 


Here the signatnre of the metric gMN is chosen to be (-I-, —, —, —, —), M is the fnndamental 
hve-dimensional energy scale, V{(j)) is a bnlk scalar held potential, Vi^ 2 ( 0 ) are brane scalar held 
potentials, g = deig^^, and g^^i, denotes the metric indnced on the branes. The space of extra 
dimension is the orbifold S^/Z 2 , which is realized as the circle of circnmference 2L with the 
points y and —y identihed. Correspondingly, the metric gMN and the scalar held 0 satisfy the 
orbifold symmetry conditions 


9^,u{x, -y) = ggu{x, y), g^^ix, -y) = -gy,i{x, y), 
gAi{x, -y) = ^ 44 ( 3 :, y), 0 (x, -y) = 0 (x, y). 


(4) 


The branes are located at the hxed points of the orbifold, ?/ = 0 and y = L, and it is assnmed 
that the SM helds with Lagrangian Lsm live on the brane aX y = L. 

For this form of the action, the vacnnm solntion for gravity and the stabilizing scalar held 
and the vacnnm solntion for the SM helds are independent. If we consider this brane-world 
model to be an indivisible theory, it is reasonable to believe that there shonld be a common 
interconnected vacnnm solntion for all these helds. To this end, let us modify action ([2]) so 
that it contains a quartic interaction of the stabilizing scalar held and of the Higgs held, which 
would connect the stabilization of the extra dimension size and the SM spontaneous symmetry 
breaking: 

S 4 >+sM = j d^x j dy (^^g^^dM<pdN(l)-V{(j))^ ^/g (5) 

- [ \/^Vi{(j))d‘^x + f ^/^{-V2{(j)) + Lsm-hp + Lint{(j),H))d^x, 

J y=0 J y=L 
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where the Lagrangian Lsm-hp is the SM Lagrangian without the Higgs potential, which is 
replaced by the interaction Lagrangian 


= , ( 6 ) 

^ being a positive dimensionless parameter. A similar interaction Lagrangian quadratic in the 
stabilizing scalar held was discussed in [T8] . 

The background solutions for the metric and the scalar held, which preserve the Poincare 
invariance in any four-dimensional subspace y = const, look like 

— dy^ = 'jMN{y)dx^dx^, (7) 

(j){x,y) = (j){y), (8) 


with denoting the hat Minkowski metric, whereas the background (vacuum) solution for 
the Higgs held is standard, 

Hvac = > ( 9 ) 


with all the other SM helds being equal to zero. 

If one substitutes this ansatz into the equations corresponding to action ([1]), one gets a 
relation between the vacuum value of the Higgs held and the value of the held 0 on the brane 
at y = L, 


0^(L) 


Mv^ 


( 10 ) 


This means that in such a scenario the Higgs held vacuum expectation value, being pro¬ 
portional to the value of the stabilizing scalar held on the TeV brane, arises dynamically as a 
result of the gravitational bulk stabilization. 

The gravitational background solution follows from the well-known system of nonlinear 
diherential equations for functions A{y),(j){y) [6], 


% + -L) = + 0" (11) 

12M\A'Y + \{V - =0 

\ {W? + ^ + V^5{y) + V25{y - L)) = -2M^ (-3A" + Q{A’f) , 


where the hrst equation is the equation for the scalar held and the other two follow from the 
Einstein equations. Here ' = d/dy. 

Suppose we have a solution A{y),(j){y) to this system for an appropriate choice of the 
parameters of the potentials such that the interbrane distance is stabilized and equal to L. 
This means that the vacuum energy of the scalar held has a minimum for this value of the 
interbrane distance. 

Now the linearized theory is obtained by representing the metric, the scalar and the Higgs 
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field in the unitary gauge as 


gMN{x, y) 

(t){x,y) 

H{x) 


'jMNiy) + 


0(l/) + 

0 

v+a(x) 

~VT~ 


V 2 IP 


huN^x, y), 


V 2 W 


f{x,y), 


( 12 ) 

(13) 

(14) 


substituting this representation into action ([T]) and keeping the terms of the second order in 
hMN, f and a. The Lagrangian of this action is the standard free Lagrangian of the SM 
(i.e., the masses of all the SM helds, except the Higgs held, are expressed in the same way 
as usual in terms of the vacuum value of the Higgs held and the coupling constants) together 
with the standard second variation Lagrangian of the stabilized Randall-Sundrum model [1] 
supplemented by an interaction term of helds / and a. The part of the Lagrangian relevant to 
the Higgs-radion mixing is 


2M3 V2 # 


1 d‘^V2 2Xv‘^^ 

^ M 


f2 , 1 2 2 

/ H-r— fa - 2Xv a 

M2 2 




(15) 


The bulk scalar held / can be expanded in KK modes. In the case under consideration we 
can, in the standard way, hnd the equations for the mass spectrum and the wave functions of 

in the formulas 


1 Cp\2 


-)■ 


1 d^V2 

2 drj>^ 


+ 


2\v‘^ 


the KK excitations of the scalar helds just by replacing 2 ~d^ ' 2 ' m 

of paper [1] (we note that in this paper the signature of the metric was (—, +, +, +, +), and for 
this reason some formulas of that paper my diher in sign from the corresponding formulas in 
the present paper). Namely, in this paper it was shown that it was most convenient to describe 
the scalar degrees of freedom of stabilized brane-world models by the held y = e~'^^^'^'^hi 4 ,{x,y) 
related to the held / by the gauge condition 


^2A 


f = -SAfi—x'. 


(16) 


The equations of motion and the boundary conditions on the branes for the held y dehne wave 
functions {yn(|/)}, corresponding to mass eigenvalues Expanding the held / in these modes, 
substituting this expansion into the second variation Lagrangian and integrating over the extra 
dimension coordinate y, we get a four-dimensional Lagrangian, in which there is an interaction 
between the modes and the Higgs held coming from the term fa. Using gauge condition 

dH, the mode decomposition of x{x,y), 


x{x,y) =^(j)n{x)Xn{y), 

n=0 


and the boundary condition for the mode wave function Xn{y) at y 
of eq. flT^ . 

(IrPV^ ^ j M 

\2 d<t''‘ 4 j' M ) - 


= L [1] rewritten in terms 
0 , 
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the field / can be expressed through Xniy) as follows: 


3M3 




\y=L — 


X I 

' A>'y=^ ~ 


CO 




n=l 


f 1<PV2 , J_ 2>^\ ,, 

2 (j)' ^ M ) ^ 


Xn{L)(^n{x), 


where we have taken into account that A{L) = 0 in order to have Galilean coordinates on the 
brane ai y = L [3]. 

Thus, the couplings of the modes to the Higgs field can be written as 


CO 

^ylan(j)n{x)a{x), 

n=l 


where we have introduced dimensionless quantities 


1 d'^V2 I 
2~d^ ' 


+ 


2Ad2^\ 

M J 


'Xn{L) — ttr, 


(17) 


which are supposed to be positive. 

These couplings are proportional to the squared masses of the modes and can be rather large 
for large n. Below we will see that it is the dimensionless quantities a„ that are of importance. 
They are proportional to the values of the wave functions of the modes on the brane and it is 
convenient to express them in terms of the ratios of the latter and ax 

Xn 
Ctn 

Xi 

These ratios must go to zero for large n in order for the cumulative effect of the radion KK 
tower to be finite. In paper |3] the scalar wave functions {Xn{y)} have been found explicitly 
in a stabilized brane-world model, where the warp factor is approximately equal to that of 
the unstabilized Randall-Sundrum model, and one can check that the quantities a„ are indeed 
positive and the ratios really fall off fairly quickly for large n in this case. It is rather difficult 
to find them explicitly in the stabilized brane-world models, where the warp factor is different 
from the exponential of a linear function, as it is in the Randall-Sundrum model. A study 
of such stabilized brane-world models has been carried out in papers Ell], and it was found 
that they may also solve the hierarchy problem of the gravitational interaction, giving rise to 
the masses of KK excitations in the TeV energy range. However, the corresponding equations 
cannot be solved exactly and should be studied numerically for each set of the fundamental 
parameters of the model, which is a very complicated and laborious task. 

Here, we will not carry out such calculations for a specific stabilized brane-world model, but 
rather give a qualitative description of the phenomena due to the interaction of the Higgs field 
with the radion and its KK tower, choosing the masses and the coupling constants in a consistent 
manner and taking into account the present-day theoretical and experimental restrictions on 
their values. 


{L) 


{L) 


(18) 
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3 Higgs-radion mixing and the effective Lagrangian 

The part of the extended SM Lagrangian containing the scalar fields linearly or quadratically 
looks like 


-'part 


- ^2Ar; ^ 

n=l n=l 


(19) 


X] K4>n(Tl; + Arp 

n=l / ^ 


n=l 


9 /If 2 /If 2 Af 2 Af 2 

+ ^Z,Z^a\ 


V 


V 




2^2 


where is the trace of the SM energy-momentum tensor and AT^ is the conformal anomaly 
of massless vector helds explicitly given by 

l^{9s)^ah 


AT,^ = 


2g^ GpG£ + 


( 20 ) 


with jS being the well-known QCD and QED [3 functions. 

The coupling of the scalar modes to the trace of the energy-momentum tensor comes from the 
standard interaction Lagrangian for the metric fluctuations [1], which results in the expression 
for the parameters in terms of the scalar mode wave functions: 


bn 


1 


fXn{L). 


2\/^' 

In the case of the lowest scalar mode, the radion, this parameter is usually denoted as 


( 21 ) 


''■“V 


( 22 ) 


where is supposed to be of the order of the fundamental energy scale M. Below we will use 
Ar as a natural scale for the interactions of the radion and for those induced by its KK tower. 
For this reason, it is also convenient to express bn in terms of A,, and an defined in (ITSD as 




(23) 


In what follows, we consider the phenomenology of the Higgs boson and the radion in the 
energy range much lower than the masses of the radion excitations. In this case, we can pass 
to a low-energy approximation for this Lagrangian by dropping the kinetic terms of the radion 
excitation helds and integrating them out, which gives the following effective Lagrangian for 
the interactions of the Higgs and radion helds with the SM helds 


Apart—ef f 


ia^tra^tr - i(2Ati^ - - i/jp? + iJ-lhh'r 

Z4„(T!; + ATp - ^a(Tl: + ATp - 

Ivj- ^ U 

9 A/f2 AL2 A/f 2 AL2 

^W-W^+a + ^ZpZ^a + + ^Z,.Z^a^ 

A* n, y 2 


(24) 
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V 


2^2- 
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where the new parameters are defined in terms of the old ones as follows: 


j2 2 2 2 

d =a^ 

n=2 


c = aiJ2^l- 

n=2 


( 25 ) 


We have already mentioned that these series should be convergent. The positive parameter 
can be of the order of because the mass term 2Au^ — should be positive, and, within 
the approach under consideration, the coupling constant of the Higgs boson self-interaction A 
can be larger than in the regular SM case. Possible restrictions on the parameters c will be 
discussed below. 

Next we turn to the mass matrix of the helds a and 0i, which looks like 

(26) 

V 2^'!®! hi / 

This matrix can be diagonalized by a rotation matrix 


Tl'^MTZ = diag{rnl,ml), TZ 


f cos 6 — sin 6\ 
ysind cosd j ’ 


where the rotation angle 6 is dehned by the relation 


tan 2d 


fil — 2Xv‘^ + d^ 


(27) 


It is not difficult to see that it is sufficient to take this angle only in the interval —vr/d < d < vr/d, 
which is one full period of the function tan 2d. 

The observable masses of the mass eigenstates are given by the expressions 


m. 


= 2\v^-d‘^ 


jji\ — 2Xv^ + d? 


'1 + 


da^ 


(/if — 2Xv‘^ + 


= 


A<;+ — 


2Xv^ + d^ 


1 + 


daf 


(/if — 2Au^ -I- d^)' 


- 1 


which are written in the form that explicitly shows that the smaller mass becomes smaller and 
the larger one becomes larger due to the mixing. The latter, in particular, means that the 
observable masses mf and cannot be exactly equal in the presence of the mixing and the 
mixing angle d is negative for mf < in accordance with the observation in paper [6]. 

The original unobservable parameters of the matrix Ad can be expressed in terms of the 
observable parameters , mf and d. The correspondence between these sets of parameters is 
one-to-one, if one puts the mixing angle d = 0 for = mf. Then the following formula for 
the parameter ai can be easily obtained: 


Gl — 


(mf — m\) sin 2d 
mf cos^ d -|- sin^ d 
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This formula, together with the second formula in ([25]), gives an expression for the parameter 
c in terms of the physical masses and the mixing angle 



(28) 


which also depends on the sum of the wave function ratios These ratios are, of course, model 
dependent, and, although they should fall off with n in order for the sum to be convergent, 
one cannot exclude that, in certain models, several first ratios can be of the order of unity. 
Thus, one can conservatively estimate this sum to be also of the order of unity0 which gives 
the estimate for the parameter c 


0 ^ C ^ Cmn.T, 


(m^ — m^) sin 26 
cos^ 9 + sin^ 9 ’ 


(29) 


which will be used in the subsequent calculations. 

The physical mass eigenstate fields h{x),r{x) are easily expressed in terms of matrix TZ and 
the original fields as 


h{x) = cos9 a{x) + sm9 (j)i{x) (30) 

r{x) = — sm9 a{x) + cos9 (f)i{x). 

The field h{x) is called the Higgs-dominated field, and the field r{x) is called the radion- 
dominated field because cos 6* > | sin6*| in the interval —7r/4 < 6 < vr/d (we recall that 6 < 0 
for < ml). The interaction of these fields with the fields of the Standard Model is given by 
the following effective Lagrangian: 


L 


h—r — 


+ 

+ 


^df.hd’^h - )mlh^ + ^d^rd^r - 

+ AT.) + + AJ..) 


— -'iljf-ifjf{cos9 h — sin^r) H-^^^“^^^^(cos^h — sin6*r) 

f V V 

— —Zfj_Z^(cos9 h — sin 6'r) -|-^lT'“lT^’''(cos6'h — sin 6'r)^ 

V ^ 

M| 


2n2 


Z^Z^{cos9 h — sin 6'r)' 


(31) 


The effective four-dimensional interaction Lagrangian (l3T]l expressed in terms of the physical 
Higgs-dominated h{x) and radion-dominated r{x) fields contains their interactions with the SM 
fields and involves only five parameters in addition to those of the SM^I: the masses of the Higgs- 
dominated and radion-dominated fields m/j and m^, the mixing angle 9, the (inverse) coupling 

^For example, in the stabilized model discussed in paper [4], the approximate wave functions of the scalar 
modes give rather rough estimates for the value of this sum in the interval (0.02,0.2) for the radion mass in the 
interval 100 GeV < rrir < 500 GeV and the mass of its first KK excitation of the order of 1 TeV. 

^We do not assume a priori with which state, the Higgs-dominated or the radion-dominated, the 125-GeV 
boson observed at the LHC is associated. 











constant of the radion to the trace of the energy-momentnm tensor of the SM helds and the 
parameter c that accommodates the contribntions of the integrated-ont heavy scalar modes. 
We would like to note here that, according to formulas fl25]) and fl2^ . the parameter c is 
proportional to its maximal value Cmax and therefore depends on the masses rrih, mr and the 
mixing angle 6, which is a particular feature of the model. 

It is interesting to examine what occurs to Lagrangian fl3T|) if the fundamental energy scale 
M goes to inhnity. In the model under consideration, it is impossible just to take the limit 
M —)■ oo because there exists relation (1T0|1 for the model parameters, which includes the vacuum 
value of the Higgs held. In the stabilized brane-world model discussed in [1], the radion mass 
was found to be proportional to 0^(L)/M^. Thus, in order to keep this mass hxed, we have 
to take the limit M —)■ oo and ^ 0 so that = const. One can see that in this case Oi, 

which is proportional to does not vanish, whereas dehned in and fl22l) goes to 

inhnity. As a result, in Lagrangian fl3ip the terms with the energy-momentum tensor that are 
proportional to l/A^ vanish, but all the other terms remain because the mixing angle must not 
be equal to zero. The situation in this case is rather similar to the one in the SM extended by 
an extra singlet scalar [191120]. However, there are diherences due to extra interactions between 
the Higgs held and the singlet scalar in these models, arising from the scalar held potentials 
that are absent in our model. 

We would also like to note here that, if we formally put the parameters Oi, c and 6 equal 
to zero, i.e., consider the case of the zero mixing, Lagrangian (ITI]) becomes just the SM Higgs 
Lagrangian plus the usual Lagrangian of the radion interaction with the SM helds via the trace 
of the SM energy-momentum tensor. In the case of a no n zero mixing, there are additional 
terms in this Lagrangian that may lead to certain changes in the collider phenomenology of the 
Higgs boson and the radion. These issues will be discussed in the next section. 

4 Phenomenological constraints from the LHC 

The ehective Lagrangian fl3ip describes the interactions of the Higgs-dominated h{x) and the 
radion-dominated r{x) helds with the Standard Model helds and with each other. A few natural 
questions arise. What model parameter regions are allowed by the present-day experimental 
data such as LEP and Tevatron searches, Higgs-like 125-GeV boson discovery and the Higgs 
signal strength measurements at various signatures at the LHC, and the limits from searches 
for the second Higgs-like boson at the LHC? Is an interpretation of the observed boson as a 
radion-dominated state is still possible for some model parameters? 

Detailed answers to such questions require a corresponding detailed phenomenological anal¬ 
ysis, which will be considered in a separate study. In this paper we present constraints coming 
from the discovery of a 125-CeV Higgs-like boson at the LHC, measurements of the signal 
strengths at various production channels and decay modes, and searches for heavier Higgs 
bosons. 

The Feynman rules needed for our study can be easily derived from Lagrangian fl3Tl) and 
are given in the table in the Appendix. We have omitted the terms coming from off-shell 
fermions since the corresponding contributions to physical processes are exactly canceled out 
by additional diagrams with contact four-point vertices coming from the hrst-order expansions 
in gauge couplings of the SM energy-momentum tensor trace as was proved in [2T] . 
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The Feynman rnles have been implemented as a “new model” in a special version of the 
CompHEP code [2^ [231123] which also inclndes a rontine for analysis of the Higgs signal 
strength in a way similar to those nsed in papers [231123] . The next-to-next-to-leading order 
(NNLO) corrections are taken into acconnt in the CompHEP compntations and in the cor¬ 
responding analysis by mnltiplying the involved vertices by correction factors for each model 
parameter point snch that the partial and the total SM Higgs decay widths and the Higgs 
production cross sections in gluon-gluon fusion (GGF) are exactly equal to those presented by 
the LHG Higgs Gross Section Working Group [271 EHj- Because of the similarity of the Higgs 
boson and the radion production and decay amplitudes including loops [2T] , the same correction 
factors have been used for the Higgs- and the radion-dominated states. 

We examine two possible scenarios, where the observed 125-GeV boson is either the Higgs- 
dominated state or the radion-dominated state. 

4.1 Higgs-dominated state at 125 GeV 

In order to understand how much room is left for the radion-dominated state, we hrst demon¬ 
strate various decay and production properties of this state. In Fig. [H the decay branching 
ratios for the radion-dominated boson are shown as functions of its mass. 

Radion-dominated state branch fr. (A,=3 TeV, c=0, sin0=O) Radion-dominated state branch fr. (A,=3 TeV, sine=0.7) 




Figure 1: Radion-dominated state branching ratios as functions of nir {jrih = 125 GeR, = 
3 TeV). The left plot corresponds to sin^ = 0 and c = 0; the right plot corresponds to 
sin^ = 0.7, and c is equal to its maximum value Cmax- 

The left plot shows the well-known branching behavior [29] of the radion without mixing 
with the Higgs boson and interacting with the SM helds only via the trace of the energy- 
momentum tensor. As is well known, the main decay mode of the light radion is the decay to 
two gluons due to the anomaly enhancement. However, once the mixing with the Higgs boson 
is included, the picture might be changed drastically because of the appearance of Higgs-like 
couplings for the radion-dominated state proportional to sin 6 and because of the influence of 
the higher scalar modes encoded in the parameter c. 

In Fig. 121 the behavior of the upper level of this parameter as a function of the radion mass 
and the mixing angle is illustrated. We recall that for 6 > 0 the radion-dominated state is 
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heavier than the Higgs-dominated state, > rrih, and vice versa ( 6 * < 0 for < rrih), and 
therefore, for 6 < 0, rUr > rrih and for 6 > 0, rrir < rrih, fhe areas shown as zero plates are not 
allowed. If one takes a rather large mixing angle, for example, sin0 = 0.7, the bb decay mode 
is dominating as shown in Fig. [T]on the right plot. Because of a compensation in the gluon- 
gluon-radion vertex between the trace anomaly part proportional to l/A^ and the mentioned 
Higgs-like part proportional to sinO/v the vertex can be very small for some particular regions 
of the parameter space. This leads to an interesting feature shown in Fig. [Uthat the decay 
branching ratio for the radion-dominated state to two gluons may be smaller than the branching 
ratio to two photons. Because of the mentioned compensations between various parts of the 


The upper level of the parameter c 



Figure 2: Two-dimensional plot of the upper level Cmax of the parameter c as a function of 
{rrir, sin 9). 

interaction vertex of the radion-dominated state with gluons, the behavior of the total width 
and the production cross sections for 77 and ZZ* modes in gluon-gluon fusion also have some 
minima, as one can see from Figs. [3] and HI From the cross sections one gets the corresponding 
signal strengths for the radion-dominated state divided by the SM Higgs cross sections. In 
our analysis we use the Higgs signal strengths for all the channels, as given recently by the 
CMS [30] and ATLAS [31] collaborations. As a result of the standard analysis, one gets 
the regions in the parameter space (the mass of the radion-dominated state and the mixing 
angle with the Higgs boson) which are still allowed. The regions allowed either by considering 
only the 77 mode, and the 77 together with the ZZ* modes are given in the left and the 
right plots of Fig. H] correspondingly, for a not-too-large mass range for the radion-dominated 
state. For both cases, all the main production processes [GGF, vector boson fusion (VBF), 
and associated production with vector bosons (VH) and with the top quarks (ttH)] are taken 
into account. The interference effects of the Higgs- and radion-dominated states are also taken 
into account, which is especially important in the case of close resonances. The contours in all 
hgures correspond to 65%, 90% and 99% best-£t conhdence level (CL) regions with less 
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than 2.10, 4.61 and 9.21, respectively, which is the standard presentation of two-parameter fits 
and was used in earlier papers (see, e.g., |2Sl 132] )• Thus, the dark shaded area corresponds to 
65% CL of the fit, the medium shaded area corresponds to 90% CL of the fit, and the light 
shaded area corresponds to 99% CL of the £t. 


Higgs-dominated state width (m^=125 GeV, Ar=3 TeV, c=c^^) 


Radion-dominated state width (m=125 GeV, Ar=3 TeV, c=c^^^) 



Figure 3: Three-dimensional plots of the total width as a function of {mr,sm6) for the LHC 
at y/s = 8 TeV and m/j = 125 GeV, = 8 TeV, c = Cmax- The left plot corresponds to the 
Higgs-dominated state h{x), and the right plot corresponds to the radion-dominated state r{x). 


99 ^ yy (LHC, Vs=8 TeV, m^=125 GeV, A^=3 TeV, c=c^^) gg ^ ZZ* (LHC, is=8 TeV, m^=125 GeV, A^=3 TeV, c=c^^) 



Figure 4: Three-dimensional plots of the partial cross section as a function of {m^, sin 6 ') for the 
LHC at ^/s = 8 TeV and ruh = 125 GeV, Ar = 8 TeV, c = Cmax- The left plot corresponds to 
the gg —?• 77 channel, and the right plot corresponds to the gg ^ ZZ* channel. 
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fit yy tag (LHC, V^=8 TeV, m^=125 GeV, A^=3 TeV, c=c^^) 


fit YY tag + ZZ* tag (LHC, \[s=8 TeV, tT^=125 GeV, A=3 TeV, 



rrir, GeV rrir, GeV 

Figure 5: Exclusion contours for the partial fit in the {rrir, sin 6) plane for the LHC at ^/s = 8 TeV and 
ruh = 125 GeH, A^. = 8TeV, c = Cmax- The dark, medium and light shaded areas correspond to 65%, 90% 
and 99% CL of the fit, respectively. The left plot corresponds to the pp —>■ 77 channel, and the right plot 
corresponds to the combined fit of the pp —^ 77 and pp —^ ZZ* channels. 


Global fit (LHC, V^=8 TeV, (t^= 1 25 GeV, A,=3 TeV, c=c^^) 



m,, GeV 

Figure 6: Exclusion contours for the combined 
fit in the (mr,sind) plane, which includes all the pro¬ 
duction processes (GGF, VBE, VH and ttH) and all 
the main decay channels ( 77 , ZZ*, lElE*, 66 , t + t “) 
for the LHG at ^/s = 8 TeV and mu = 125 GeF, 
Kr = 3 TeV, c = Cmax- The dark, medium and light 
shaded areas correspond to 65%, 90% and 99% CL 
of the fit, respectively. The white cross marks the 
best-fit point. 


Global fit h-resonance (LHC,t[s=8 TeV, rr^=125 GeV, A,=3 TeV, c=c,„^) 



m,, GeV 

Figure 7: Exclusion contours in the heavy ra- 
dion mass parameter for the combined fit in the 
{rrir, sin 0 ) plane that comes from the signal strengths 
at 125 GeV for all the production processes (GGF, 
VBE, VH and ttH) and all the main decay chan¬ 
nels ( 77 , ZZ*, IVIV*, 66 , T+r“) for the LHC at ^/s = 
8TeV and ruh = 125GeF, Ar = 3TeF, c = Cmax- 
The dark, medium and light shaded areas correspond 
to 65%, 90% and 99% CL of the fit, respectively. 
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In Fig. [5l one can see that the 77 mode alone allows the presence of the radion-dominated 
state for masses both below and above the 125 GeV. The ZZ* mode gives some additional 
restrictions on the parameter space especially in the mass region closer to the Z-boson pair 
threshold, where the cross section is increased. When all the modes are taken into account, the 
allowed parameter space region is reduced further, which is shown in Fig. [ 6 l 

For larger values of the radion mass, the influence of this parameter on 125-GeV signal 
strength is very small, resulting in the allowed area shown in Fig. [71 However, in the case of a 
large resonance mass region, one has to take into account the exclusion limits given very recently 
by the GMS [33] and ATLAS [3l| collaborations, coming from searches for heavy Higgs bosons. 
If the radion-dominated state has a large mass, it would contribute as a resonance in this region. 


Global fit r-resonance (LHC,Vi=8 TeV, n^=125 GeV, A,=3 TeV, c=c^^) 

O 
C 

0.6 

0.4 


0.2 

0 

- 0.2 

- 0.4 

- 0.6 

200 300 400 500 600 700 800 900 1000 

m,, GeV 

Figure 8 : Exclusion contours in the heavy radion 
mass parameter for the combined ht of the CMS 
and ATLAS exclusion regions at a high mass range 
by the resonance production and decay of the radion- 
dominated state with corresponding heavy mass in the 
{rriri sin plane. The dark, medium and light shaded 
areas correspond to 65%, 90% and 99% CL of the fit, 
respectively. 



Global fit (LHC, (s=8 TeV, m =125 GeV, A,=3 TeV, c=c^^^) 
c 

to 0.6 

0.4 

0.2 

0 

- 0.2 

- 0.4 

- 0.6 

200 300 400 500 600 700 800 900 1000 

m,, GeV 

Figure 9: The region in the (?7ir,sin6*) still allowed 
by the common fit of the CMS and ATLAS exclu¬ 
sion limits for heavy Higgs searches and restrictions 
from the fit of the influence of those parameters on 
the signal strengths at the point 125 GeV. The dark, 
medium and light shaded areas correspond to 65%, 
90% and 99% CL of the fit, respectively. 



From the comparison of the radion cross section computation and the experimental limits by 
GMS and ATLAS by performing the corresponding £t, one gets a rather large allowed area in 
the {rrir, sin 6 *) parameter plane, as shown in Fig. |H1 When taken together, fits of the data from 
both the signal strengths and from the direct searches for heavy Higgs resonances significantly 
restrict the allowed region for the heavy radion mass to a small area at very high masses, as 
demonstrated in Fig. [9l 

Obviously, if one considers larger values of the parameter A^, the cross section of the radion- 
dominated state gets smaller and the allowed region for such a state is increased. This is 
demonstrated in Fig. [TOl where the parameter A^ is chosen to be 5 GeV with all the other 
parameters being the same as for the considered 3-TeV case. 
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Global fit (LHC, Vs=8 TeV, rT|^=125 QeV, A,=5 TeV, o=c,^^) 


Global fit (LHC, is=a TeV, m^=125 GeV, A,=3 TeV, c=c^^) 



trir, GeV 

Figure 10: Exclusion contours for the combined 
fit in the (nir, sin 0) plane that includes all the pro¬ 
duction processes (GGF, VBF, VH and ttH) and all 
the main decay channels {jj, ZZ*,WW*,bb,T~^T~) 
for the LHC at ^/s = 8 TeV and mu = 125 GeF, 
Kr = 5 TeV, c = Cmax- The dark, medium and light 
shaded areas correspond to 65%, 90% and 99% CL 
of the fit, respectively. The white cross marks the 
best-ht point. 



GeV 

Figure 11: Exclusion contours for the combined 
fit in the (ruh , sin 0) plane that comes from the sig¬ 
nal strengths at 125 GeV for all the production pro¬ 
cesses (GGF, VBF, VH and ttH) and all the main de¬ 
cay channels ( 77 , ZZ*, WW*, bb, t+t“) for the LHC 
at ^/s = 8TeV and nir = 125 GeF, = 3TeF, 
c = Cmax- The dark, medium and light shaded ar¬ 
eas correspond 65%, 90% and 99% to CL of the fit, 
respectively. The white cross marks the best-fit point. 


4.2 Radion-dominated state at 125 GeV 

Now let us consider briefly the case where the radion-dominated state has mass 125 GeV, and, 
using the same analysis strategy. And the regions for the mass of the Higgs-dominated state 
and the mixing angle allowed by two signal strengths. For this case, we have carried out all the 
corresponding calculations and have drawn the plots similar to those in Figs. |2H3 which we 
will not present here for the sake of brevity. We only dwell upon the result of the analysis 
that is presented in Fig. [11] 

As one can see, such a scenario is strongly disfavored in comparison with the previous case. 
The mass of the Higgs-dominated state might still be very close to the mass of the radion- 
dominated state in a wide range of mixing angle. Such a possibility exists in both scenarios. 

5 Conclusion 

In the present paper we have considered the Higgs-radion mixing arising in stabilized brane- 
world models due to merging the mechanism of stabilization of the extra dimension size and 
the Higgs mechanism of spontaneous symmetry breaking on the TeV brane. We have discussed 
phenomenological restrictions on model parameters coming from coupling measurements and 
searches for heavy bosons at the LHC. This mixing is, of course, similar in many aspects to the 
one arising due to the Higgs-curvature term on the brane. However, an important difference is 
the presence of an extra coupling at low energies of the Higgs-dominated held to the trace of 
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the SM energy-momentum tensor originating from the coupling of this held to the heavy scalar 
states of the radion KK tower. 

In order to study the physical consequences of the Higgs-radion mixing in stabilized brane- 
world models, we derived the effective Lagrangian and gave a qualitative description of the 
phenomena, taking consistent values for the masses, the coupling constants, and the mixing 
angle, which are useful for comparing the results of our calculations with the experimental data 
at the LHC. It turned out that, though the interaction of an individual higher excited scalar 
state with the Higgs held may be weak, their cumulative effect on the Higgs-radion mixing 
may be observable. If they give noticeable contributions to the parameter c, it leads to certain 
changes in the collider phenomenology of the Higgs boson. A similar contribution of the directly 
unobservable higher tensor KK modes to scattering processes was discussed in [35] . 

In the framework of this model, we have studied two a priori possible scenarios: The scalar 
state discovered at 125 GeV is either a Higgs-dominated state or a radion-dominated state. 
In our analysis, we used measurements of the 125-GeV Higgs signal strengths by the ATLAS 
and GMS collaborations and exclusion limits obtained in searches for heavy Higgs-like states. 
Our results show that the interpretation of the 125-GeV scalar state as a Higgs-dominated 
state is the preferred one, although the radion component in this state can be rather large. 
Depending on the value of the radion coupling constant A^, the allowed regions for the mass 
of the radion-dominated state have been found. It turns out that the radion-dominated state 
can either have a mass close to 125 GeV or a mass close to the TeV range. The allowed regions 
somewhat increase with the growth of the radion coupling constant A^. 

We have also shown that the interpretation of the 125-GeV scalar state as a radion- 
dominated state is not completely excluded by the two leading signal strength measurements, 
though in this case the restrictions on the allowed masses of the Higgs-dominated state are very 
stringent. The mass of the Higgs-dominated state can be close to 125 GeV, which is in accor¬ 
dance with our analysis of this state at 125 GeV. Thus, in the considered model the presence 
of two nearly degenerate states close to 125 GeV is a very probable scenario. 
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7 Appendix: Feynman rules 


In Table 1, the constants and the functions of the momenta are explicitly given by the following 
expressions: 5 qcd = 7, 62 = Y’ ~ “V’ ~ sin^ -|- c • cos^, Cr = cos 9 — c ■ sin 6 ', 

Fw = -(2 + 3yw + ^ywi‘2 - yw)f{yw)), Ft = yt{l + (1 - yt)f{yt)), yi = 4m,V(2pi • P 2 ), 


f{z) 


sm 


log 


-1 

1 + 

1 - 
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z > 1 
z <1. 
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